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ABSTRACT 


We propose a mechanism for the nonthermal ion heating 
observed in the solar wind at lA.U. based on an electrostatic, 
short wavelength. Instability between the ions in the observed 
colliding plasma streams. The modes lying on a plane perpen- 
dicular to the magnetic field thermalize most of the differential 
energy. The model suggests local electrostatic turbulence, that 
a-particles are heated more than protons (T^ ^Tp) , a maximum 
proton temperature Tp s? lo\ and for the bulk speed U possibly 
a Tp = aU + b relationship. These predictions are consistent 


with observations. 



I. 


INTRODUCTION 


Direct satellite observations during the last 9 years have 
led to a remarkably detailed and reasona>>ly complete knowledge of the 
interplanetary plasma (e.g; Strong et al 1966; Neugebauer and Snyder 
1966J Hundhausen 1968; Burlaga and Ogilvie 1969; Parker 1969) . 

Equilibrium hydrodynamic models of the coronal expansion while" fairly 
successful in predicting the ovserved mass flow properties underestimate 
the proton temperatures in the active wind by more than one order of 
magnitude (Sturrock and Hartle (1966)) . In addition/ the distribution of 
the internal energy among the various ion species indicates clearly 
the existence of non-thermal heating mechanisms. Proper identification 
of these mechanisms seems to be the key in understanding the observed 
thermal properties of the interplanetary plasma. 

Sturrock and Hartle (1966) first suggested that the required 
energy for ion heating may be contained either in magnetohydrodynamic 
waves originating near the Sun and propagating outward through the solar 
wind plasma or- in the differential motion of solar wind streams. In 
the latter case the Kelvin-Helmholtz Instability is the mechanism by which 

the differential energy is converted into long wavelength turbulence. 

* •'/ 

Following this Coleman (1968) argued that regardless of the source of 
the waves, it is likely that they cascade down into shorter wavelengths 
by non-linear processes, much as in hydrodynamic turbulence, till they 
transfer their energy to the ions by ion cyclotron damping. Jokipii and - 
Davis (1969) argue that if the waves were of solar origin they would damp 
out before reaching the earth's orbit as shown by Barnes (1967). Sub- 
sequently assuming that such waves are generated by the colliding plasma 



streams, suggest that they heat the plasma directly by resonant inter- 
actions and most of the wave energy does not cascade to smaller eddiss, 
as proposed by Coleman. 

The colliding stream structure of the wind, was discovered 
by Snyder et al (1963) and is at times associated with the magnetic sector 
structure of Ness and Wilcox (1966). Sarabhai (1963) noted that as solar 


rotation brings a sector of higher wind velocity beneath one of lower 
velocity, the fast stream overtakes the slower one. J okipii and Davis 
(1969) argue that it is the collision between such streams that generates 
the waves, rather than the shear between them (Coleman 1968). However 
in either view the observed low frequency and large wave-length fluctu- 
ations are associated with the "anomalous" ion heating, which should 
therefore be a large scale phenomenon (Parker 1969) . 

The recent observational analysis of the Explorer 34 data by 


Burldga and. Ogilvie (i969) seems to contradict various aspects of the 
proposed above non-thermal ion heating mechanisms. 

Lack of correlation between negative velocity gradients (as 
viewed from the earth) with the obse'rved non-thermal ion heating seems 
to rule out the hypothesis of the Kelvin-Helmholz instability. 

The observed turbulence, as opposed to fluctuations, occurs 


only in- "patches" and is well correlated with local proton heating. This 
is consistent with the interpretation of data from Mariner 2 by Neugebauer 
and Snyder (1966) and from the Vela by Hundhausen et al (1967) , who placed 


an upper limit oh’ the 'wavelength of" the observed turbulence of the order 
5-50 km (Hundhausen, 1968). Such a limit implies that the magnetic field 
does not participate in the turbulent or collective motions, which should 
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therefore be electrostatic in nature. 

Finally strong correlation of local proton heating with large 
positive velocity gradients, is consistent with the idea that heat is 
generated by the collision of fast streams with slower moving plasma. 

These observations indicate that non- thermal ion heating is 
a localized rather than large scale phenomenon, electrostatic in nature, 

.. > r»- ^ 

consistent with the hypothesis of electrostatic instabilities in col- 
liding streams and not connected directly with the large scale hydro- 
magnetic turbulence. 

It is the purpose of the present note to propose an ion heating 
model consistent with the above picture. The thermalization mechanism 
is a short wavelength electrostatic instability between the counter- 
streaming ions in the colliding streams. Under some local conditions 
specified in the 'next section modes with wave vectors (k) on a plane 
perpendicular to the magnetic field are unstable, and non-linear theory 
indicates that they result in thermalizing the differential energy per- 
pendicular to the direction of the magnetic field. An extensive account 
of the linear and non-linear theory of the instability supported by 
computer experiments has been published recently by Papadopoulos et al 
(1971) and Wagner et al (1971). Below, following a brief summary of 
the instability theory, we present arguments which suggest that the 
proposed model is consistent with the various observational characte- 
ristics. 

II. THEORETICAL BACKROUND • 

Several studies of ion instabilities between plasma streams ap- 
peared lately in the literature due to their importance in high Mach 
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number shocks and controlled thermonuclear research (Papadopoulos 

et al 1971; Wagner et al 1971; Landau 1971). We present here a short 
review of the results pertaining to unstable modes on a plane perpen- 
dicular to a magnetic field, which appear to be important in the 
problem under consideration. For a detailed account we refer the inte- 
rested reader to Papadopoulos et al (1971) and Wagneiret al (1971). 

The system shown in fig. la is linearly unstable if 

u $ 2.5 A + 3 V. (1) 

• r A 

where is the local Alfven spead (V^ = B/[4nnM]^) and 3 the ratio 

2 

of the electron to the magnetic nressure (3 = nk„ T /B /8tt) . The growth 

B e 

ratfe and the wavenumber for the most unstable electrostatic mode is 

2 

- CJ 

... = 1/2 03^ k^ = 3/4 ~ (2) 

mom 2 

u - - 

r 

2 2 2 2 

where m = o)./2[l + co /f2 ] , m (to.) the electron (ion) plasma frequency 
o 1 e e e 1 

and the electron cyclotron frequency. The range of the electrosta- 

tically unstable wavenumbers k is given by 

'2(0^ 

0 < k < \ (3) 

u 

r 

2 2 2 

subject to the condition k c /to^ (1 + 3) 5 1 • Electromagnetic instabili- 
ties do not result in ion thermalization and can therefore be neglected 
for the present case. 

Non-linear theory shows that the instability results in thermalizing 
the relative ion streaming energy producing a single-humped ion distribution 
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function moving with the velocity of the center of mass of the system 
(fig. lb) and temperature 


T. . 

^ -3 • (4) 

1/2 Mu . 

r 

The stabilization mechanism is ion trapping in the large amplitude 
electrostatic ion turbulence which occurs at a flucttiation energy level 
c„ such ‘ 


-2 
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These thoretical results have been confirmed by computer simulation 

4 

using a code which follows the orbits of 2 x 10 electrons and ions. Fig. 

1 (a,b) shows the evolution of the ion distribution function and ion 
phase space for a typical numerical experiment, with parameters relevant 
to the solar wind. 

In the discussion presented above we assumed equidensity ion beam 
streaming perpendicular to the direction of the magnetic field B . One 
can generalize the results for cases of oblique streaming by replacing the 
relative streaming speed u^ by its projection u^ sin 0 in the direc- 
tion perpendicular to the magnetic field. This can relax the instability 

\ 

condition so that streams with relative velocity u $ 2V. /l + B/sin 0 

X k 

can be unstable. However the temperatures achieved will be lower than 

2 

givenby eq. (4) by a factor sin 0 . Therefore streamers along the mag- 
netic field lines will not be heated. The effects of unequal stream den- 
sities were discussed in Papadopoulos et al (1971) . The essential conclusion 
is that while for density ratios close to unity the ion temperatures will be 
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of the order of magnitude given by eq. (4), For large density ratios only 
the low density stream will be therraalized and will thus form a hot ion 
tail in the main ion distribution function. In all cases the average ion 
heating is related to the increase in the bulk velocity (Au) over that 
of the slow plasma with a relation of the form 

= aAu (6) 

where a will be a function of the density ratio and the angle 9 . 

' In closing this section we mention that the instability has been 
recently identified as producing ion heating in two laboratory experiments 
(Davis et al 1971; Dean et al 1971). 

Solar Wind Ion Heating 

We examine next whether the implications of the proposed model 
are consistent with the observations of the non-thermal properties in the 
solar wind. We'^oncentrate on observations related to the turbulence and 
overall ion heating, rather than the anisotropic features such as the ratios 
Tj| /T_j^ of the parallel to perpendicular temperatures which should be accounted 
for in terms of internal instabilities of electromagnetic nature (Parker 
1963; Scarf et al 1967) and are independent of the electrostatic mechanism 
producing local ion heating. These observations, as reported in an excellent 
review article by Hundhausen (1968), and a later paper by Burlaga and Ogilvie 
(1969), are: 

(1) Proton temperatures : The observed proton temperatures in the active 

wind are of the order of 3 x 10^ with an upper limit of 9 x 10^ °K 
(Hundhausen 1968). These are more than one order of magnitude higher than 



4 

the 10 °K expected from the two fluid model (Sturrock and Hartle 1966 ) . 

The proposed mechanism can explain both the observed temperatures and the 
upper limit by independent considerations. As indicated by eq. (4) in order 
to achieve temperatures of the order of 3 ^ 10^ °K one needs relative speeds 
between the streamers of the order of 140 km/sec, which definitely falls in 
the range of relative speeds expected and observed. ’^The maximum temperature 
of the ions due to the proposed thermalization mechanism can be found by. com- 
byning eqs. (1) and (4) and is. given by (T ) 1/4 6.25 (1 + 3) 1/2 V^. 

^ 1H3.X if A 

Such temperatures of course would appear only if u = 2.-5/l + 3 V. which 

r A 

is the largest relative velocity for the system to be unstable. Since for 
the active wind at lAU 3 .5 and V, 'u 80 km/sec. we find that (T„) 

10^ ®K . It is remarkable that his upper limit set by the values of magnetic 
field and the density as measured in the active wind rather than by speculations 
about the differential energy of the stream, is in such a good agreement with 
the proton temperature observations. 

(2) Temperature distribution among the various ion species; It has been 

observed that while in the quiet wind the temperatures of the various ion 
species present are almost the same, in the active wind they are greatly 
different and tend to be proportional to their masses (e.g., the helium 
to proton temperature is 4) indicating equality of their thermal 

velocities. This seems to be one of the strongest arguments in favor of a 
two stream electrostatic instability. Prom eq. (4) we find that the stabi- 
lization temperature ' T. of the ion species j with mass M. follows the law 
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Notice that the constant depends only on the relative streaming , 

which is presumably the same for all species. The observed distribution 
of Internal energy among the various ion species is an absolutely natural 
consequence of the proposed mechanism. . 

(3) Scale and wavelengths of the turbulence: Measurements of flux vs. 
angle or energy per change peaks show a substantial spread whic-h might be 
the result of turbulent or collective motions. Such an interpretation of the 
data as discussed by Neugebauer and Snyder (1966) places an upper limit of 

50 km on the wavelength of the' turbulent motions as seen in a reference frame 

moving with-^the-plasma flow. Similar considerations on the data transmitted 

by the Vela satellite lower the wavelength limit to 5 km (Hundhausen et al 

1967). This is consistent with the Explorer 34 observations (Burlaga and 

Ogilvie 1969) which indicate that turbulence occurs only in "patches" 

and should be distinguished from fluctuations which are a large scale effect. 

As was discussed in the previous section the electrostatically unstable 

2 2 2 

wavenumbers are such that k c /o) i 1 . Since in the active wind at lA.U 
• e 

c/o)^ 1 km, we would expect the maximum wavelength to be X 'v 2Tr/k 6 km, 

which seems to be consistent with the observations. In the proposed model 

we would expect the .electrostatic ionic turbulence to be imbedded within, the.. 

large scale hydromagnetic fluctuations produced either by the colliding 

streams, in regions where they are electrostatically stable (u > 2.5/1 + .8V. ) 

r A 

or by internal instabilities due to anisotropies in temperature. 

(4) Flow. speed to proton temperature correlation: High bulk velocities 

in the solar wind- have been- found to correlate well with large proton tern- - 

peratures (Neugebauer and Snyder 1966). In -particular Burlaga and Ogilvie . 

% 

find an empirical relation of the form T^ = aU + b to be valid over a wide 
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variety of solar wind conditions. The Sturrock-Hartke model predicts 

quite accurately the observed Tp-U relation for the quiet wind. Thus 

if this remarkable relation is actually a fundamental property of the solar 

wind flow, the increase in temperature (ATp 'v Tp) should be accompanied 

% 

by an increase in the bulk speed AU of a functional form TpUAU . This 

also is a consequence of the proposed mechanism since as one can see from 

. the theory of the instability part of the relative streaming energy goes to 

random energy and part to increase the bulk speed over that of the slower 

stream. As explained in Papadopoulos et al (1971) the proportionality constant 
% • 

a in the Tp = aAu relationship is a function of the density ratio of the 
interacting streams and the angle 0 . Although the exact numerical Tp-U 
relationship cannot be found from our simple model, we believe that the pro- 
perty of the instability to increase the bulk speed proportionally to the 
thermal speed is a 'positive indication. 

(5) Ion tails; We mention finally that tails in the proton distribution 
function such as expected when low density streamers collide with high density 
streamers, are observed regularly in the wind (Hundhausen 1968). In addition 
the fact that electrons and protons are heated by different mechanisms (Bur- 
laga and Ogilvie 1969) is consistent with the proposed model. 

As mentioned at the beginning of the section besides the properties 
discussed above, the proton temperature exhibits an anisotropy Tp^ /Tpj_ 2 . 

We feel that this is not connected with the electrostatic hon-thermal heating 
mechanism which occurs over such short distances that the direction of the _ 
magnetic field over the interaction times depends more on the local fluctu- 
ations than in equilibrium field considerations. This is more so since in 
regions of colliding streams magnetic flux conservation requires a temporary 
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and local enhancement of the field in the direction perpendicular to 
the streams. Thus if the ratio Tp|j/Tpjj_ is controlled by internal 
electromagnetic instabilities (Parker 1963; Scarf et al 1967) would 

therefore expect such instabilities which thend to restore isotropy, to 
follow the local proton heating, producing long wavelength fluctuations, - 
such as observed. 

Before closing this section we mention that-' one can ..apply similar 
considerations to the thermalization of the differential motions described 
by Coleman (1968) which are due to large scale fluctuations in the magnetic 
field and the streaming velocity. This can provide an ion heat source even 
in the quiet wind. 

Summary and discussion; 

In the foregoing the hypothesis that the non-thermal proton 
heating and turbulence in the solar wind Is due to a short wavelength, 
electrostatic instability between the counterstreaming ions in colliding 
plasma streams has been tested. It has been shown that the dominant non- 
thermal properties of the solar wind are consistent with the predictions 
of the proposed mechanism. The picture emerging is that of localized 
proton heating and short wavelength electrostatic turbulence embedded within 
large scale magnetohydrodynamic fluctuations. It was pointed out that 
following the proton heating, instability related to the emerging tempe- 
rature anisotropies might create long wavelength fluctuations. Additional 
proton heating might arise from thermalizing differential motions due to 
the large scale fluctuations. 
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Fig. la 



Initial ion phase space and distribution 
function 

Fig. lb Ion phase space distribution function at 
thermalization. 
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